Prospects for constraining the shape of non-Gaussianity with the scale-dependent bias 
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We consider whether the non-Gaussian scale-dependent halo bias can be used not only to constrain 
the local form of non-Gaussianity but also to distinguish among different shapes. In particular, we 
ask whether it can constrain the behavior of the primordial three-point function in the squeezed limit 
where one of the momenta is much smaller than the other two. This is potentially interesting since 
the observation of a three-point function with a squeezed limit that does not go like the local nor 
equilateral templates would be a signal of non-trivial dynamics during inflation. To this end we use 
the quasi-single fleld inflation model of Chen & Wang [JQ] as a representative two-parameter model, 
where one parameter governs the amplitude of non-Gaussianity and the other the shape. We also 
perform a model-independent analysis by parametrizing the scale-dependent bias as a power-law on 
large scales, where the power is to be constrained from observations. We find that proposed large- 
scale structure surveys (with characteristics similar to the dark energy task force stage IV surveys) 
have the potential to distinguish among the squeezed limit behavior of different bispectrum shapes 
for a wide range of fiducial model parameters. Thus the halo bias can help discriminate between 
different models of infiation. 

I. INTRODUCTION 

The study of the deviation from Gaussianity of the initial conditions set by inflation is one of the most active fields 
of research in cosmology today. This is for a good reason since it is a potential observational handle on the interactions 
of the inflaton. Indeed, a free field is Gaussian and all its correlation functions are fixed by Wick's theorem in terms 
of its two-point correlation function 

(C(^i)C(fc2)) = {27TfS{h + k2)Pc{k) , (1) 

where C is the comoving curvature perturbation, and = 27r^A^(fc/fcp)^"=^^)(l/fc^). The spectral index and 

amplitude of the scalar power spectrum have been measured by the WMAP satellite Q to be Ug = 0.968 ± 0.012 
and A^ = (2.430 ± 0.091) x 10~^ at the 68% confidence level, and kp is an arbitrary pivot scale there chosen to be 
kp = 0.002 Mpc^^. The interactions of the field induce non-trivial higher-order correlation functions which if observed 
can serve to discriminate between different models of inflation. In the context of cosmology, Non-Gaussianity refers to 
the study of the generation of these higher-order correlation functions in different models of the early universe, their 
modification through the non-linear evolution of the later universe, and how observations may be used to measure 
them. Of these, the three-point function is expected to be the most accessible to observations; it can be written as 

(C(fci)C(fc2)C(fc3)) = {2TTfS{k, + k2 + k-s)Bc{ki, fc2, fca) , (2) 

where is called the bispectrum. Early work focused on the study of the three-point function given by the following 
phenomenological prescription 

c = Q + lfhLic'-{0), (3) 

where Q is a Gaussian random field. This generates a bispectrum with the following shape 

4°^(fci,fc2,fc3) = l.f^r[Pdki)Pc{k2)+cycl.] . (4) 



2 



This type of bispectrum is called local since it stems from a local redefinition of the field. Note that it is divergent 
in the limit in which one of the momenta is much smaller than the other two, often called the squeezed limit; this 

divergence goes as (C(fci)C(fe)C('Z)) l/?^- 

The three-point function generated by the simplest single-field slow-roll models with a standard kinetic term was 
first computed by 0, [l] , who found it to be too small to be accessible to observations in the foreseeable future. A 
"large" non-Gaussianity can be generated by introducing new ingredients in the model such as additional light fields 
(see e.g. [l]), a non-trivial kinetic term (see e.g. 0)j or features in the power spectrum (see e.g. [1]). Each of these 
models generates a characteristic type of non-Gaussianity. Furthermore, within the context of single-field slow-roll 
inflation, the observation of the three-point function could potentially fix the coefficients of the leading terms in the 
effective Lagrangian for the perturbations of the inflaton [§, [T^] . 

The best constraints on the primordial three-point function to date come from the WMAP satellite observing the 
cosmic microwave background radiation (CMB) Q. A non-trivial three-point function has yet to be observed, but 
the Planck satellite promises to greatly improve on current bounds Another promising probe comes from the 

observation of large scale structure in the universe (LSS). Competitive constraints on the non-Gaussianity have been 
found [III using the fact that a non-Gaussianity of the local type, Eq. induces a characteristic scale, rcdshift and 
mass dependence on the halo bias [isl - fisj . 

In this paper we wish to study the power of LSS surveys in the near future to constrain the shape of non-Gaussianity 
through the observation of the halo bias. The halo bias is sensitive to a configuration in which one of the three Fourier 
modes is much smaller than the other two, the so-called squeezed limit, where the long mode is the scale of observation. 
The local non-Gaussianity has a divergent squeezed limit and gives a contribution to the large-scale halo bias that 
goes like ^ In section HIl we summarize the information about infiation contained in the squeezed limit of the 

bispectrum. It has been shown 0, that the bispectrum of single-field infiationary models is suppressed in the 

squeezed limit by two powers of the small momentum, thus giving a negligible contribution to the halo bias unless 
there are features in the power spectrum close to the relevant scales. Models with several light fields may generate a 
bispectrum which has a squeezed limit that behaves like the local non-Gaussianity thus giving potentially interesting 
effects to the halo bias in future surveys, as studied for example in [l9| . In order for the three-point function to have 
a non-standard behavior in the squeezed limit, some non-trivial dynamics must be present during inflation. One such 
model is the "quasi single-field" inflation model of [l|, [Ij , which is the one we use as an example of the power of a 
future survey to constrain the behavior of the three-point function in the squeezed limit. It consists of a light inflaton 
coupled to an isocurvaton with a mass which is close to the Hubble scale during inflation. In section Hill we summarize 
how the presence of primordial non-Gaussianity affects the halo bias in the peak background split formalism. In 
section HV] we describe the survey setup we consider, with the characteristics of the specific Dark Energy Task Force 
stage IV survey. In particular, we use a survey similar to that used in (20j and we refer the reader to the survey 
description given on that reference. In section |V] we present our main result, namely the forecast of the constraints 
that such a survey would be able to impose on the parameters of the quasi-single field inflation model, which can also 
be thought as a parametrization of the squeezed limit of a three-point function. We also perform the same analysis 
for a power- law parametrization of the non-Gaussian bias on large scales. Finally, in section IVTl we draw conclusions 
and discuss the results. 



II. MODEL OF NON-GAUSSIANITY 



As we will review in the following section, the non-Gaussian modiflcation of the halo bias is mainly sensitive to the 
squeezed limit of the bispectrum in which one of the momenta, here denoted by q, is much smaller than the other 
two. It is thus of interest to consider what information can be gained about inflation from constraints on the squeezed 
limit of the three-point function. 

Generic single-fleld models of inflation must satisfy the following consistency relation 

(C(fcl)C(fc2)C(g)> -(2^)3^(fci + k2+ q)P{ki)P{q) 
This relation was first derived in Refs. while Ref. [l^ showed that the corrections are suppressed by the square 



-(n.-l)+0^ 



(5) 
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of the small momentum (i.e., the second term m the square brackets). This implies that all single field models of 
inflation generate a bispectrum which in the squeezed limit goes like (C'^) ^ 1/q (^). This was understood in reference 
[Tsj to be a consequence of the fact that, under the assumption of adiabaticity, correlation functions of C are endowed 
with a SO{A, 1) symmetry which is non-linearly realized by C- 

The exact momentum dependence of the bispectra as predicted by different models of inflation is too cumbersome 
to be used in CMB data analysis. What is often done is to look for a bispectrum shape, often called a template, which 
has a high overlap with the exact bispectrum and use it to analyze the data. Two such templates which are often used 
are the equilateral one [2l| and the orthogonal one [l0|- They are similar to shapes generated by different operators 
in the effective Lagrangian for single field inflation 0,113] • However, one must be careful when using these templates 
with the scale-dependent halo bias to put constraints on the models which inspired them. Indeed, the orthogonal 
template has a squeezed limit going like {(^) ~ which is in contradiction with the above-mentioned fact that all 
single field models of inflation must have a squeezed limit going like {(^^) ~ 1/q (^). One can still view studies using 
the orthogonal template and the halo bias as examples of the usefulness of this method in constraining intermediate 
shapes. Nevertheless, at the time of writing, no model in the literature generates a bispectrum which has both a large 
overlap with and a squeezed limit like that of the orthogonal template. 

For models with several light fields which can contribute to the generation of primordial perturbations none of the 
above arguments hold, and the bispectrum can generically be large in the squeezed limit. For example, curvaton 
models (see e.g. 0]) predict a bispectrum which is very close to having a shape like that of the local template, and 
which goes like (C^) l/?^ in the squeezed limit. This is one of the reasons why observing a large non-Gaussianity 
of the local type is compelling, as it would rule out all single-field models of inflation. However, there is still more 
information contained in the behavior of the bispectrum in the squeezed limit. Reference [l7j argues that in multi- 
field models, when all the fields are much lighter than the Hubble scale during inflation and are thus scale invariant, 
the squeezed limit of the bispectrum is proportional to the power spectrum of the long mode, i.e., cx with 
corrections which are again quadratic in the small momentum. Thus, in single-field models of inflation and in models 
with multiple light fields the bispectrum always goes as or l/g in the squeezed limit. 

If one wishes to obtain a non-standard behavior in the squeezed limit, the alternative is to consider a multi- field 
model where at least one of the fields has a mass close to H, which is high enough to invalidate the arguments of 
the previous paragraph, but low enough so that it can't be trivially integrated out. An interesting realization of this 
is the model of Ref. [H, 0: termed "quasi-single field inflation". It consists of two fields, one being a light infiaton 
and the other an isocurvaton with a mass of order H. The infiaton and isocurvaton are coupled through a turning 
trajectory in field space, which allows for the isocurvaton to have a sizeable contribution to the three-point function of 

Reference [l| computes the bispectrum for such a model, and suggest the following template as a good approximate 
description 



where is the Neumann function, a is fixed to be 8, and v = ■\/9/4 — m? / H'^, with m being the mass of the 
heavy isocurvaton. The factors in front fix the normalization of the template to be equal to the local template 
at the equilateral point fci = ^2 = ^3- The amplitude /nl is related to the third derivative of the isocurvaton 
potential, and the angular velocity in the turning trajectory as /nl = oi{y)P^ ^^"^ (Oq/ HY{—V"' / H), where a{v) is 
the function plotted in figure 8 of Ref. [l| . 

In the squeezed limit, the template of Eq. ^ goes like (C'^) ^ l/g'^^^^'^; thus different values of v give different 
behaviors in the squeezed limit. For v = 1.5, the template behaves like the local one in the squeezed limit; for v = 0.5 
it behaves like the orthogonal template in the squeezed limit; other values of v give intermediate behaviors. One 
therefore expects future surveys to constrain the value of v using the halo bias. 



^ In order to prove this, one must assume that the time variation of C, outside of the horizon is suppressed at least as . This is the case 

for most models of inflation present in the literature. 
^ It is worth noting however that Reference [To| gives in their Appendix a template for the orthogonal shape which does have the correct 

behavior in the squeezed limit. 
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The halo bispectrum is actuaUy not sensitive to the exact squeezed limit in which one of the Fourier modes is 
infinitely much smaller than the other two. This means that even if all single field models must satisfy the consistency 
relation of equation ([S]) in the exact squeezed limit, there might be single field models with features in the bispectrum 
close to the actual scales relevant for observations such that they can still produce a large effect in the halo bias. Some 
work in this direction has been done in Ref. [j^ . Another possibility is that the bispectrum breaks scale invariance, 
again in this case the model must satisfy the consistency relation of equation ([5]) but can still have an effect on the 
halo bias, see Ref. [2^. Refs. [H, HH] very recently argued that a modified vacuum state can have a feature near to 
the squeezed limit that induces a scale dependence of the halo bias which is stronger than that generated by the local 
template. 

In this paper we wish to forecast the constraints that a future survey with the characteristics of the Dark Energy 
Task Force stage IV [2^ can put simultaneously on v and /nl as defined in equation In the context of the 
quasi-single field model, a constraint on v translates into a constraint on the mass of the heavy isocurvaton. More 
generally, one can see this as an estimate of the power of future surveys to constrain the behavior of the squeezed 
limit of the three-point function. In section Hill we also introduce a parametrization of the non-gaussian modification 
to the halo bias on large scales in order to perform this analysis in a model independent way. All of this is potentially 
interesting since if a bispectrum with a squeezed limit that does not go like l/g or l/g'^ is observed, it would signal 
the presence of heavy fields or other non-trivial dynamics during inflation. 

Let us now estimate the values of /nl that are allowed by current observations. Though no dedicated data analysis 
has been performed for the template of Eq. (jS]), it has a sizable overlap with the local template. To be more 
quantitative, define the "cosine" between two bispectra as 

cos(Bi,B2) = -=£y^= where ■ = V Si(fci, fca, ^3)52(^1, fca, fcg) . (7) 
V Bi ■ B1B2 • i>2 I I J, 

One can then use current constraints on the local, equilateral and orthogonal templates Q to estimate the values of 
/nl compatible with current data analyses: 

/NL«/NrC0s(S<;,S;„e) + ... , (8) 

where the dots indicate other templates which have been analyzed; we don't use them here since the tightest constraints 
come from those on the local template. Indeed, the cosine between the template under study, Eq. ([5]), and the local 
template, Eq. is always > 0.6. We plot the estimated l-cr constraints in figure [TJ Note that this is only a rough 
estimate, since the true constraints coming from the CMB data would require a dedicated data analysis, which is 
beyond the scope of this article. One can also improve on this estimate by using a spherical cosine rather than the 
flat one defined in Eq. ([7]); we expect the results to be of the same order of magnitude as those plotted in figure [T] 



III. MODELING OF THE NON-GAUSSIAN HALO BIAS 



In this section we wish to summarize the computation of the modification of the halo power spectrum in the presence 
of non-Gaussian initial conditions. Let us begin by introducing some notation. 

The comoving gauge curvature perturbation in terms of which the initial conditions set by infiation are often 
expressed, is related to the linear dark matter density contrast 5 at redshift z through the transfer function T and 
the linear growth factor D normalized to one at redshift zero 

5{k) = k^T{k)D{z)C{k) , (9) 

where we have written this equation in Fourier space. These factors contain the evolution of the gravitational 
potentials once the perturbation enters the horizon^. The transfer function as written here can be readily computed 



^ Note that a diflfcrent definition of the transfer function is often used, such that it is normalized to one on large scales. We do not 
normalize to one, but rather use equation l[9j as the definition, which is the one that appears directly in the output file of CAMB at 
redshift zero (though one must be careful with the factors of h). 
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FIG. 1: Estimated WMAP 7 constraints on /nl at l~a coming from the overlap between the quasi-single field inflation template 
of Eq. Q and the local template of Eq. (jU as defined in Eq. ([7|. 

using publicly available numerical codes like CAMB [27t . 

In the high-peak formalism, one considers a halo to form when a region of space has an average density contrast 
larger than some threshold Sc- Let us define this average as 



(2^ 



W{kR)S{k) , 



(10) 



where ly is a window function that averages out modes with a wavelength larger than the region under consideration. 
We will take this filter to be a top-hat in real space, which in Fourier space takes the form 



W{kR) 



{kRy 



: (sin kR — kR cos kR) . 



(11) 



The critical density contrast threshold dc is often estimated by computing the time at which a perfectly spherical halo 
completely collapses and extrapolating the linear density to that point. In an Einstein-de Sitter space one obtains 
Sc ~ 1.686, while in ACDM this estimate receives a small correction and is slightly redshift dependent. The scale 
R over which one averages is related to the mass of the halo formed by M = Air/S^lmPcR^ where pc is the critical 
density of the universe pc = STJq/SttG. Equations ((9)) and pO)) suggest the following definition 



Muik) = k^T{k)W{kR) 



(12) 



where R is expressed in terms of M as above. Thus for example the variance of the filtered density contrast can be 
computed as 



'M 



d^k 

(27r)3 



■Mlj{k)P^{k) 



In the presence of non-Gaussian initial conditions, the (Eulerian) halo power spectrum P{k) is modified 

P(fc) = (fe(f'+A6)2p,(fc), 



(13) 



(14) 
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where g denotes the Gaussian case. It was reahzed in Refs. |12| - |15| that this correction to the power spectrum can 
have a characteristic scale dependence that makes it relatively easy to observe. One can estimate this correction with 
an arbitrary bispectrum shape using the peak background split formalism as done originally in Ref. [l^ 

M, = f + '-^] . (15) 

MM[k) y D(z) dlncTM J 

where the second term in the parenthesis was computed in Refs. (28l. [29|. and is due mainly to the modulation of the 
variance with mass due to the non-Gaussianity. The "form factor" T for an arbitrary bispectrum is defined as 

F{k,M)^ ^ 2 „ I dkiklMMiki) f dfiMMiJk"^ + + 2kikfi)B^{k, ki,Jk^ + kf + 2fcifc/^) . (16) 



In reference [28| equations ([T5|) and (|16p were found to agree well with simulations [30| including a local, orthogonal 
and equilateral templates for the bispectrum of the initial conditions. This agreement has been seen to be good up 



to a scale oi k ^ 0.1 /iMpc [31 



In the integral of equation (jl6p . the scale k corresponds to the scale of observation which for the survey we will 
study varies between k ^ 0.003 /iMpc^^ and** k ~ 0.1 /iMpc^^. On the other hand, for a halo mass of order 
M = 10"^^ Mq/i"^, the filter W{kiR) suppresses the integral at scales fci > 3 /iMpc~^, while the transfer function 
and the factor of k\ in the integral suppress the integral at scales fci ^ 0.1 /iMpc~^. The measure of the integral 
will thus peak at fci ^ 0{1) Mpc^^. Therefore the integral will receive a sizable contribution only from the squeezed 
configuration of the momenta appearing as the arguments of the bispectrum in which k is at least an order of 
magnitude smaller than fci. We will use equations (jlSp and (jl6p to compute the non-Gaussian modification to the 
halo bias induced by a template of the form given by equation (|6|). 

At very large scales the non-Gaussian modification to the halo bias probes very squeezed configurations. In such 
a regime Ab{k,M) often behaves as a power law in fc as can be seen from equations ^T5\i and ((T6| . Indeed one 

expects that for a bispectrum going as (C^) l/q", the bias behaves as Ab{k, M) ^ l/fc°~^. In order to perform a 
model-independent (though perhaps rough) analysis of the power of future surveys to constrain the behavior of the 
halo bias at large scales, we use the following parametrization (as suggested by |25| ) 

A6(fc,Af) = /SL^, (17) 

where A{M) is a mass-dependent normalization defined such that the bias is the same as that in the presence of a local 
non-Gaussianity at a scale of fc = 0.03 /iMpc"^^, and /3 is a parameter to be constrained. Note that /nl of equation 
(pl and 9''''^ normalized differently, so that in order to compare the results obtained when using the template of 
Eq. ([6]) and those obtained using the parametrization (flT)) one should multiply f^^ roughly by a factor of 8. 



IV. SETUP 



We wish to estimate the capability of a future survey satisfying the requirements of the Dark Energy Taskforce 
stage IV |26| . to constrain the shape of the bispectrum. For this we assume a survey similar to the one studied in 
reference [20] ■ We present the characteristics we assume for such a survey in table [H We also use the results of Ref. 
[20| for the average number of galaxies available to the survey, the average mass of the galaxies observed, and the 
Gaussian bias; let us summarize them here. 



* In previous analyses that did not consider the correction to the halo bias given by the second term in the right hand side of Eq. III6I I. 
the maximum k considered was 0.03/i/Mpc. This more conservative cutoff ensured that the adopted bias description was correct (on 
these scales the correction is negligible) but reduced the total signal-to-noise. 
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Sky coverage 2 x 10 square degrees 

Minimum redshift 0.5 
Maximum redshift 2.1 
Typical galaxy halo mass 10^^ Mq/i-i 

TABLE I: Description of the Dark Energy Task Force stage IV survey used. 

For the Gaussian bias wc use the results reported in figure 4 of Ref. [l^l, who in turn use the results of Ref. 
[3^ . There the bias is estimated using the halo model, a semi-analytic model for galaxy formation, and assuming a 
spectroscopic selection based on Ha emission with a threshold given for a survey similar to Euclid. 

The modification of the halo power spectrum due to non-Gaussianity will depend on the halo mass. Here we take 
the typical masses of the haloes of the galaxies observed to be 10^^ Mq/i"^ as in [l^l, see their figure 3. We compute 
the non-Gaussian modification to the halo bias evaluated at this fixed mass. We don't expect the errors induced by 
fixing the mass in equation (|16|) to be of qualitative importance for our results. 

The number of galaxies at a given redshift can be computed from the halo mass function n{z^ M) (the number of 
haloes in a differential mass interval at a given redshift), and the first moment of the halo occupation distribution 
{Ng\M) (very roughly the probability that there be Ng galaxies in a halo of mass M) 

/>oo 

ng{z)= / n{M,z){Ng\M)dM , (18) 

JMg 

where Mg is the minimum halo mass for a galaxy observed by the survey for which we again follow Ref. (20l | and take 
Mg = 10^^ Mq/i"^, see their section 4.4. Since the average number of galaxies is important only in computing the 
shot noise, which we have verified has a small effect, we neglect the non-Gaussian correction to the mass function and 
use the Sheth and Tormen mass function l33i . As for the first moment of the halo occupation distribution, we follow 
Ref. and use the following expression [SJ, 

{Ng\M) = Ng,,,{^—^ , (19) 

and the parameters 9, Mq and Ng Q depend on the type of galaxy considered. We assume here typical galaxies to be 
blue, for which A'^^o = 0.7 and = ii M < Mq and 9 = 0.8 otherwise and take Mg = 4 x 10^^ MqH''^ as in section 
4.1 of Ref. [2^. Let us stress that the average number of galaxies affects only the shot noise which we verify to be 
small and our results will not be very sensitive to the assumptions made in computing it. 

V. FORECAST 

Wc work under the approximation that the Likelihood for the (Eulerian) halo power spectrum is a Gaussian centered 
around a "fiducial" model for which the values of the parameters of interest /nl and ly are fixed at /nl and D 

= (20) 

z, ^ p 

where AP = P — P\ j^^.p is the deviation from the fiducial model. The standard procedure is to assume the behavior 
of the halo power spectrum to be nearly linear on the parameters of interest (here called generically 0i), so that 
the Likelihood function is a Gaussian also in terms of those parameters. One can then estimate the variance and 
covariance of a future survey by computing the Fisher information matrix on the parameters 
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This is a good approximation when one is interested in a region in parameter space around the fiducial model which 
is small with respect to the typical variation of the power spectrum with the parameters, i.e., when the resulting 
estimated uncertainties are small. However, in the model we are studying the variation of the power spectrum with 
respect to v is highly non-linear and this approximation need not hold. We will instead compute the A^^ of each 
point in parameter space 



where we sum over the different redshift shells and Fourier modes available to the survey. This quantity is in general 
costly to compute when dealing with a high-dimensional parameter space. However, we will be interested only in 
varying /nl and v since our main focus is the measurement of the non-Gaussianity and its behavior in the squeezed 
limit. This implies that we keep other cosmological parameters fixed at the WMAP 7 best fit value s 3 , and the 
actual uncertainties on /nl and ly are expected to be somewhat larger than our estimates. Reference [36[ estimated 
this increase in the errors to be mild in the case of a local non-Gaussianity (around 10% ^ 30% depending on the 
survey if the number of relativistic neutrino species is kept fixed). This can be understood as being due to the fact 
that no other cosmological parameter induces a scale dependence and a redshift dependence on the halo bias as the 
one induced by a local non-Gaussianity. The same argument holds for the range of u we will study, and we thus expect 
the increase in the estimated uncertainty on /nl and v due to the uncertainty on other cosmological parameters to 
be small. 

For an infinitesimal shell in Fourier space with a width Ak and effective volume Veff the relative error can be 
estimated by counting the number of Fourier modes 



2 



P J 47rfc2AfcT4///(27r)3 " 
The effective volume is the volume corrected by taking into account the shot noise 



(23) 



where ng is the average number of galaxies, V is the volume of each redshift shell, and we've explicitly stated the 
dependence on redshift and Fourier mode. Equation (|22p can now be written explicitly 



0.1 h/Mpc since it is the scale up to which we trust the calculation 
of the non-Gaussian modification to the bias [3l|, and the minimum k is computed as kmin ~ 1-k lV{z,Yl^. Wc divide 
the redshift range in shells with a width Az = 0.1. The halo power spectrum and average number of galaxies are 
given respectively in equations <\1A\ and (fT8|) . The deviation of the halo power spectrum from the fiducial model was 
computed by taking the difference between the halo power spectra of the given model and the fiducial one computed 
using equation ([T5|) to model the effect of non-Gaussianity. 

In figure [2]we show the Ax^ < 2.3 regions corresponding to 68.3% joint confidence levels for different values of the 

fiducial parameters v and /nl- For i> = 0.5, which would correspond to a squeezed limit going like (C^) "^^^ l/?^, the 
constraints on v will be broad, and even the detection of non-Gaussianity would be challenging, requiring very large 
values of /nl in mild tension with current CMB constraints. As D increases, the effect on the bias becomes stronger 
and the survey becomes more sensitive to both /nl and ly. For D = 1.5, corresponding to a shape that behaves like the 
local template in the squeezed limit, the survey becomes sensitive to /nl = O(l0), and is able to put constrains on 
i>. Recall that we plot the 68% 2-parameter joint confidence regions, a similar confidence interval for a 1-parameter 
only model -i.e., if ly were to be kept fixed- would correspond to Ax^ < 1; thus to compare ours with other analyses 
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where v is kept fixed, one can "slice" the contours at the preferred v value, but should interpret the resulting error as 
similar to a 90% confidence region. Note that for large values of v and /nl the regions resemble ellipses, signaling the 
fact that the errors are small relative to the variation of v and the Fisher approximation holds. For smaller values of 
V and /nl this is not the case and the confidence regions adopt more complex shapes. 
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FIG. 2: Regions in the /nl, v parameter space defined in Eq. (|6| satisfying Ax^ < 2.3, corresponding to the 68.3% confidence 
level. In the squeezed limit this model behaves as {Q^)q^o = ^/q'^'^^^^- We show such regions for several fiducial models, 
showing that the uncertainties decrease as the fiducial value of /nl becomes larger and the fiducial value of v approaches 1.5, 
which corresponds to a bispectrum shape that behaves like the local template in the squeezed limit. 

In figure [31 we present the same analysis performed for the model-independent parametrization of equation (jl7p . 
Here we restrict ourselves only to large scales where this parametrization can be relevant, and thus integrate up to 
kmax = 0.03 /i Mpc~^. The amplitude /nl of the template of Eq. ([5]) and /^^ defined in Eq. P7)) are normalized 
differently, and one should multiply /^^ by a factor ~ 8 in order to compare the results. Thus, the values of /^^ 
in figure [3] arc taken to be small for ease of comparison with figure [2] We present the results for (3 = 1,2, and 3 
corres pon ding respectively to the behavior of the orthogonal template^ the local template, and the models studied in 
Rcfs. [24l. [25t . Again the confidence regions tend to ellipses for larger values of the amplitude /^^ ^-iid power /?, and 
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the results are consistent with those obtained for the quasi-single field inflation model in figured] 




2.6 2.8 3.0 3.2 3.4 



FIG. 3: Forecasts for the parametrization A6 = f^j^A/k^. Green regions correspond to points in the f^-^, P parameter space 
satisfying Ax^ < 2.3, corresponding to the 68.3% confidence level. Note that /^j^ is defined such that the non-Gaussian 
modification to the halo bias, eq. (|15|l . is the same as in the presence of a local non-Gaussianity at a scale k = 0.03 ft Mpc~^ , 
thus/NL(i^= 1.5) ~8/^l(/? = 2). 



VI. SUMMARY AND DISCUSSION 

We studied the possibility that a future galaxy survey with the characteristics of the Dark Energy Task Force stage 
IV might be able to constrain the shape of the primordial non-Gaussianity through observations of the halo bias. 
A non-vanishing primordial three-point function induces a modification of the halo bias that is sensitive mainly to 
the squeezed limit in which one of the three momenta is smaller than the other two. Single field models produce a 
bispectrum with a squeezed limit going like (C^) ^ 1/<Z, where q is the small momentum, while models of inflation 
with multiple light flelds can produce a large non-Gaussianity of the local type, i.e., going like (C^) ~ ^/q^ in the 
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squeezed limit. Thus, the observation of a squeezed hmit with a behavior which is different to these two possibihties 
might signal the presence of some non-trivial dynamics during inflation. Such can happen for example in the "quasi- 
single held" inflation model, which consists of a light inflaton and a massive isocurvaton coupled through a turning 
trajectory. It can produce a shape going as 

(C(fci)c(fc2)c(gi) (26) 

where v = yj^j^ — m? /H'^. 

We estimate the ability of a survey such as that described in section IIVI to constrain the squeezed limit of the 
primordial three-point function. For this purpose we perform forcasts on the uncertainties on a measurement of the 
amplitude and shape of the non-Gaussianity generated by the quasi-single held inflation model (see figure [5]) and of 
a phenomenological model-independent parametrization of the non-Gaussian halo bias at large scales A6 = f^-^A/k^ 
(see figure 131). 

Since the variation of the halo power spectrum with the power v ov (5 is large for small values of the amplitue /nl 
or J the Likelihood function in terms of these parameters is not expected to be a Gaussian. We cannot then use 
the Fisher information matrix approach, and instead compute Ax^ at each point in parameter space. We report our 
findings in figures [5] and [31 We find that the forecasted uncertainties depend quite strongly on the fiducial model. 
If the fiducial model is taken to be similar to a local non-Gaussianity {D = 1.5 or /3 = 2), the survey will be able 
to put tight constraints on the parameters of the model, even for /nl = 0(10). As P or /3 decreases the forecasted 
uncertainties become larger, requiring larger values of /nl or /^^ fo^' non-Gaussianity to be observable. This is to be 
expected since a milder divergence in the squeezed limit translates into a milder scale dependence in the non-Gaussian 
modification to the halo bias. 

Though the most stringent limits on the primordial non-Gaussianity to date come from CMB observations, the 
CMB has a limited capacity in constraining the shape of the three-point function. Indeed, it is not sensitive to the 
physically interesting squeezed limit, and would not distinguish between two shapes which have a large cosine but 
with different squeezed limits. In this sense halo bias observations are complementary to CMB observations and might 
open a new window to the physics of infiation. It would still be of interest however to find the constraints on the 
amplitude of the quasi-single field infiation template coming from CMB data. 

Another potentially interesting probe of the squeezed limit is the possibility of observing the /i distortion of the 
CMB spectrum recently proposed in Ref. j37| . A primordial three-point function induces a correlation between the 
fi distortion and the temperature which is sensitive to the very squeezed limit. It would be interesting to study the 
power of such a probe to constrain /nl and v. 

Finally, let us comment on ways in which one can increase the signal to noise ratio in our analysis. If the analysis 
is repeated for a photometric survey like LSST, which can observe galaxies with fainter magnitudes and thus higher 
redshift, one expects the gain in volume to reduce the errors by a factor ^ 2. Refs. (ssl - lioj argue that with the use of 
several differently biased tracers and an optimal weighting it is possible to reduce the errors on a local /nl by at least 
an order of magnitude. We expect a similar reduction to take place in our case, greatly improving the sensitivity of the 
survey. Another way to reduce the uncertainties is that proposed by Ref. (4l| who argue in favor of cross-correlating 
pixels in order to extract more information from the survey; they find that for a local non-Gaussianity the reduction 
of the uncertainties is again of an order of magnitude. If the uncertainty on /nl can be brought to be of 0(1), general 
rclativistic effects on the halo bias become relevant and might then be observable (e.g., [12] and refs therein). 

While this paper was being completed, we became aware of a similar analysis being carried out by another group, 
see Ref. [i^. We have carried out a rough comparison of our results and find them to be consistent. We are grateful 
to the authors for coordinating with us the publication on the ArXiv website. 
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